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Abstract 

Light strongly interacting super symmetric particles may be treated as partonic constituents of nu- 



(N 



^ ■ cleons in high energy scattering processes. We construct parton distribution functions for protons 

m ■ 

Tij- . in which a light gluino is included along with standard model quark, antiquark, and gluon con- 

^O \ stituents. A global analysis is performed of a large set of data from deep-inelastic lepton scattering, 

O 

^T ' massive lepton pair and vector boson production, and hadron jet production at large values of trans- 

r~| ' verse momentum. Constraints are obtained on the allowed range of gluino mass as a function of 

' , the value of the strong coupling strength as{Mz) determined at the scale of the Z boson mass. We 

(— I \ find that gluino masses as small as 10 GeV are admissible provided that as{Mz) > 0.12. Current 

hadron scattering data are insensitive to the presence of gluinos heavier than ~ 100 — 150 GeV. 
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I. INTRODUCTION 

Relatively light strongly-interacting fundamental particles may be considered as con- 
stituents of nucleons. The nature of these constituents and their experimental effects become 
evident when the parent hadrons are probed at sufficiently short distances or, equivalently, 
sufficiently large four-momentum transfer Q. The charm quark q = c and the bottom quark 
h are treated appropriately as constituents of hadrons in situations in which Q > rriq, where 
rriq is the mass of the heavy quark. Other strongly interacting fundamental particles may 
exist, as yet undiscovered experimentally, with masses lying somewhere between the bottom- 
and top-quark masses. One example is a relatively light gluino: a color-octet fermion and 
the super symmetric partner of the massless spin-1 gluon. For our purposes, we define a 
"light" particle to have a mass less than 100 GeV. In this paper, we explore the effects that 
a color-octet fermion would have on the parton distribution functions of nucleons, with a 
view toward establishing whether the set of hard-scattering data used in global analysis may 



already place significant constraints on the existence and allowed masses o 



QQ 



such states, 
as a concrete 



In our investigation, we use a light gluino from supersymmetry (SUSY) 
example, but our analysis and conclusions should apply as well to the case of a color-octet 
fermion of whatever origin. As constituents of hadrons, these color-octet fermions share the 
momentum of the parent hadron with their standard model quark, antiquark, and gluon 
partners. The distribution of light-cone momentum fraction x carried by constituents is 
specified by parton distribution functions (PDFs) as functions of both x and the scale Q of 
the short-distance hard scattering. The process-independent PDFs are essential ingredients 
for obtaining normalized predictions of rates for hard-scattering reactions at high energies. 
A simultaneous analysis of a large body of scattering data (global analysis) provides strong 
constraints on the magnitude and x dependence of the PDFs. 

In perturbative quantum chromodynamics (QCD), the existence of a color-octet fermion 
and its couplings to the standard model constituents alter the coupled set of evolution 
equations that governs the functional change of the parton distributions as momentum is 
varied. Gluinos have different renormalization group properties from those of the quark and 
gluon constituents, and the contributions of fermions in the color-octet representation are 
enhanced strongly. Within the context of broken supersymmetry, squarks (scalar partners of 
quarks) may also be relatively light, particularly those of the third generation, the bottom- 



and top-squarks [3|, y, |5|, |6[. In the study reported here, we include a gluino in our PDF 
analysis, but we neglect possible contributions from other hypothesized supersymmetric 
states with masses above a few GeV, such as bottom squarks. As explained in Sec. II, the 
effects of squark contributions on the current data are much less important than those of 
gluinos. The approximation of retaining only the light gluino contribution simplifies the 
calculations while retaining most of the relevant physics. 

In a global analysis of hadronic data, a large sample of data is studied (about 2000 points) 
from a variety of experiments at different momentum scales. The data set included in our 
study is the same as in the recent CTEQ6 [2| study done within the context of the standard 
model. The data come from deep-inelastic lepton scattering, massive lepton pair and gauge 
boson production, and hadron jet production at large values of transverse momentum. We 
apply the methodology of the next-to-leading order (NLO) CTEQ6 analysis to explore the 
compatibihty of a light gluino with the large set of hadronic data. Methods developed 
recently for the analysis of uncertainties of PDFs jsl Isl, lnj, luJ, ll2, ll3, llJ, ll^ allow us to 
obtain quantitative bounds on the existence and masses of glui nos from a global analysis. 
In early PDF analyses within the context of light gluinos [la, ll2l llS| , a gluino with a mass 5 
GeV or less was found to be consistent with the data available at that time. A more recent 
study jl9| disfavors a gluino with a mass 1.6 GeV or less. The much larger sample of the 
data in the modern fit and improved understanding of PDF uncertainties make it possible 
to derive more precise bounds. 

Light superpartners influence the evolution with scale Q of the strong coupling strength 
as{Q)- The constraints we obtain on the gluino mass from a global analysis depend signifi- 
cantly on the value of the strong coupling strength as{Mz) that is an ingredient in the global 
analysis. In Sec. II, we begin with a brief review of the dominant experimental constraints 
on Us and consider the changes that may arise if supersymmetric particles and processes 
are present. Further discussion of experimental constraints on as{Mz) may be found in the 
Appendix. We describe in Sec. II. B how we implement the NLO evolution of the PDFs, 
while including the gluino degree of freedom at leading order (LO). Once supersymmetric 
particles are admitted, they contribute to hard scattering processes either as incident par- 
tons and/or as produced particles. We therefore specify the hard scattering matrix elements 
that describe supersymmetric contributions to the rate for jet production at large transverse 
momentum. In Sec. Ill, we present and discuss the results of our global fits. The discrimi- 



nating power of our analysis depends crucially on the inclusion of the Tevatron jet data in 
the fit. The inclusion of a light gluino in the PDFs removes momentum from the gluon PDF 
at large x, tending to depress the contribution from SM processes to the jet rate at large Ex- 
However, the effect is compensated partially by a larger value of as{Mz), slower evolution 
of as that makes as^Ex > Mz) larger than in the standard model, and by contributions to 
the jet rate from production of SUSY particles in the final state. 

Our conclusions are summarized in Sec. IV. We find that the hadron scattering data 
provide significant constraints on the existence of gluinos whose mass is less than the weak 
scale ~ 100 GeV. A large region of gluino parameter space is excluded by the global analysis 
independently of direct searches or other indirect methods. The quantitative lower bounds 
we obtain on the gluino mass must be stated in terms of the assumed value of the the strong 
coupling strength as{Mz). For the standard model world-average value as{Mz) = 0.118, 
gluinos lighter than 12 GeV are disfavored. However, the lower bound on rrig is relaxed to 
less than 10 GeV if as{Mz) is increased above 0.120. 

II. Us, PARTON DENSITIES, AND HARD-SCATTERING SUBPROCESSES 

The presence of a light gluino g and/or a light squark q modifies the PDF global analysis 
in three ways. First, the gluino and squark change the evolution of the strong coupling 
strength as{Q) as the scale Q is varied. Second, the gluino and squark provide additional 
partonic degrees of freedom that share in the nucleon's momentum and affect the PDFs 
of the standard model partons, e.g., via the channels g ^ gg and q -^ qg. Third, gluino 
and squark contributions play a role in the hard-scattering matrix elements for the physical 
processes for which data are analyzed and fitted. We discuss each of these modifications in 
the following three subsections. 

A. Modified evolution and values of as{Q) 

The expansion of the evolution equation for as{Q) as a power series in as{Q) is 



When supersymmetric particles are included, the first two coefficients in Eq. (^ are (see, 
e.g., Ref. Q) 

2 1 

Po = 11 - 3^/ - 2ng - -nj, (2) 

and 

38 11 13 

fJi = 102 - —Uf - ASn-g - —nj + —ngUj, (3) 

where Uf is the number of quark flavors, rig is the number of gluinos, and n? is the number 
of squark flavors. Equation (J2|) shows that, to the leading order, one generation of gluinos 
g contributes the equivalent of 3 quark flavors to the QCD /9-function. The effect of one 
squark flavor is equivalent to one- fourth of the contribution of a quark flavor. In our work, 
we henceforth neglect the possibility of a light squark contribution to the /^-function and 
limit ourselves to the effects of a light gluino. Inclusion of a light bottom squark changes 
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I22I. The modifled coefficients 



the running of a^ slightly, compatible with current data 

/3o and /3i for rig = 1 and n? = are implemented in our numerical calculation to full NLO 

accuracy. 

In our global flt of hadron scattering data, the allowed range of the gluino mass rrig 
depends strongly on the assumed value of the strong coupling as{Mz). Therefore, it is 
important to understand the current experimental constraints on as{Mz) from sources other 
than hadron scattering data. A combined analysis of all Z-pole data within the context of 
the standard model, carried out by a working group of members of the four CERN Large 
Electron- Positron Collider (LEP) experiments and the SLAG SLD experiment [23], obtains 
the value as{Mz) = 0.1187 ± 0.0027.^ This value is but a shade greater than the often- 
quoted standard model world-average value as{Mz) = 0.1183 ± 0.0027 [24] obtained from 
a variety of determinations of as{Q) at different momentum scales Q. On the other hand, 
the value of as{Q) extracted from data in the context of supersymmetric contributions 
can be different from the value obtained in standard model flts. Some of the assumptions 
made in a standard model analysis are modifled by the presence of the supersymmetric 



contributions J6, 
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and the value of as{Q) at Q = Mz obtained in standard model 



analyses may have to be be revised. A recent estimate J27[ of SUSY effects on directly 
1 See Table 16.2 of Ref. 0. 



measured as{Mz) finds values in the interval (0.118 — 0.126) ± 0.005, where the variation in 
the central value arises from uncertainty in the magnitude of SUSY-QCD corrections to the 
Z-boson decay width from processes such as Z ^ bbg and Z -^ bbg. Further discussion of 
the evolution of as{Q) in the presence of light supersymmetric states may be found in the 
Appendix. 

In a general purpose CTEQ fit, as{Mz) is fixed to be equal to its world-average value, 
determined from a combination of the r-lepton decay rate, LEP Z pole observables, and 
other measurements. As discussed above, this value may change in the presence of light 
superpartners. To explore fully the range of strong coupling strengths compatible with 
the global fit, we perform a series of fits in which as{Mz) is varied over a wide range 
0.110 < as{Mz) < 0.150. We then determine the values of rrig and as{Mz) preferred by the 
global fits. 

B. Implementation of a gluino in the NLO evolution of parton distributions 

In the construction of parton distributions, we include a light gluino and omit squark 
contributions. A squark enters parton splitting functions only in combination with another 
rare particle, and these splittings are characterized by smaller color factors than in the gluino 
case. We incorporate the gluino sector into the PDF evolution package used to build the 
CTEQ6 unpolarized parton distributions jlf. 

The standard procedure for extracting parton distribution functions from global QCD 
analysis is to parametrize the distributions at a fixed small momentum scale Qq. The 
distributions at all higher Q are determined from these by the Dokshitzer-Gribov-Lipatov- 
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Altarelli-Parisi (DGLAP) evolution equations 28|, |29|, |30| . The agreement with experiment is 
measured by an effective x^? which can be defined by x^ = Sexpts Xn^ ^^ by generalizations of 
that formula to include published systematic error correlations. The PDF shape parameters 
at Qq are chosen to minimize x^ ^ind obtain the "best fit" PDFs. 

We choose the starting value Qq for the QCD evolution equal to the smaller of the gluino 
mass rrig or charm quark mass rric- At the scale Q = Qq, the only non-perturbative input 
distributions are those of the gluons g and light {u, d, s) quarks. Non-zero PDFs of the 
gluinos and heavy quarks (c, b) are generated radiatively above their corresponding mass 
thresholds. In the CTEQ6 analysis, Qq coincides with the charm quark mass: Qo = ""^c = 1.3 



GeV. Therefore, for rrig > rric the input scale Qo = 1.3 GeV is the same as in the CTEQ6 
study. We use the CTEQ6 functional forms for the input PDFs of the standard model 
partons at Q = Qq, but the starting values of the parameters are varied in order to obtain 
acceptable fits to the full set of scattering data. 

The prescription for Qq allows us to investigate the possibility of gluinos lighter than 
charm quarks {rrig < rric). We include fits for gluino masses 0.7 < nig < 1.3 GeV by choosing 
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rrin 



Such super-light gluinos may be generated both via perturbative and nonper- 



turbative mechanisms, and, in principle, an independent phenomenological parametrization 

must be introduced for the gluino PDF to describe nonperturbative contributions. Our 

prescription for the region rrig < rric provides a particular model for such an input gluino 

mrametrization, similar in its spirit to the dynamical parton distributions of the GRV group 

as well as the procedure used in earlier light gluino analyses |l8[ and 19|. ^ 

In the presence of a light gluino, the DGLAP equations must be extended to account for 



the new processes. The coupled evolution equations take the form 
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i=u,d,s,... 



(5) 



Here S(x,(5^), g{x,Q^), and g{x,Q'^) are the singlet quark, gluon, and gluino distributions, 
respectively; qi^XjQ"^) and qi{x,Q^) are the quark and antiquark distributions for flavor i. 
The splitting functions Pij{x) may be found in the literature J32l |. 

The inclusion of a gluino in the evolution equations complicates the calculation substan- 
tially. To achieve acceptable accuracy, evolution of the quarks and gluons must certainly be 
done at next-to-leading order accuracy. However, without a substantial loss in accuracy, we 



^ In Ref. 19f|, the PDFs are obtained from the procedure described here, but at LO and without inclus 
of the jet production data, for nig = 0.5 and 1.6 GeV. 
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can simplify the overall calculation by evaluating the gluino contributions to leading order 
accuracy only. We use the following prescription: 

1. Evolve the ordinary quarks and gluons at NLO (so that the splitting functions Pss, 
-Psg, -Pgs, and Pgg are evaluated to order 0(a^)). 

2. Evolve the gluinos at LO (so that the splitting functions Pgg, Psg, Pgg, Pg^,, and Pgg 
are evaluated to 0{as))- In particular, at LO (and in the absence of the squarks), 

3. For the evolution of as, use the full NLO (C(a^)) expression, including the effect of 
the gluino. 

In this prescription, the evolution is fully accurate to NLO except for the gluino splitting 
kernels. Were we interested in a process dominated by gluino contributions, we might need 
a NLO representation of the gluino PDF, g{x, Q"^). However, the impact of the gluino PDF 
is minimal for the inclusive data in the global analysis, since g{x, Q"^) is much smaller than 
the quark and gluon PDFs (cf. Figs. ^ and E))- As a result, the gluino plays only an indirect 
role. Its presence modifies the fit in two ways: 

1. The gluino alters as{Q), thereby modifying the evolution of ordinary quark and gluon 
PDFs. 

2. The gluino carries a finite fraction of the hadron's momentum, thereby decreasing the 
momentum fraction available to the gluons and standard model quarks. 

Regarding item (1), we compute the effects of the gluino correctly by using the exact 
NLO beta function that includes SUSY effects. Therefore, the only shortcoming of our 
prescription is with respect to item (2). We describe correctly the NLO mixing between the 
quarks and the gluons, but the less consequential mixing of the standard model partons and 
the gluino is correct only to leading order. In the energy range of our interest, the gluinos 
carry a small fraction (^ 5%) of the proton's momentum. The neglected NLO corrections 
to this small quantity are further suppressed by a factor of ag/vr. They are comparable in 
magnitude to the NNLO corrections for the standard model splittings, which are suppressed 
by a'^/n'^. The uncertainty introduced by the omission of the NLO gluino splittings is 



comparable to that due to the NNLO corrections for other particles, and it may be ignored 
in the present NLO analysis. 

Figures Q and |21 illustrate and support our assumptions. They show the momentum 
distributions for the gluons, singlet quarks, and gluinos, respectively. In these figures, the 
strong coupling strength as{Mz) is equal to 0.118 (the value of as{Mz) assumed in the 
CTEQ6 analysis). The momentum scale Q is 10 GeV in Fig. ^ and 100 GeV in Fig. |21 The 
abscissa (x-axis) is scaled as x^^^ in these plots of the dependence on the momentum fraction 
X. The distributions are obtained from the light gluino (LG) fits for nig = 10 GeV, with or 
without the inclusion of the Tevatron jet data. For xg{x,Q'^) and xS(x, Q^), we also show 
the corresponding distributions from the best-fit set CTEQ6M of the CTEQ6 analysis. 

With Q = nig = 10 GeV, the gluino density xg{x,Q^) = 0. Nevertheless, the effects 
of gluino contributions on the fit at Q > 10 GeV force a change in the gluon and quark 
distribution functions from their standard model values at Q = nig = 10 GeV, as shown in 
Fig. E Once Q is evolved to 100 GeV, Fig. |2l shows a nonzero momentum distribution for 
the gluinos and a persistent change of the gluon and quark densities from their CTEQ6M 
standard model values. 

The figures demonstrate two important features. First, the magnitude of the gluino 
distribution is much smaller than the gluon and quark distributions. This large difference 
justifies the assumptions that contributions are small from scattering subprocesses with 
initial-state gluinos, and that NLO gluino contributions may be omitted in our analysis. 

Second, the presence of the gluino depletes the gluon distribution at x ^ 0.05. The effect 
on the singlet distribution is less pronounced. The gluinos take their momentum (3.7% of the 
proton's momentum at Q = 100 GeV for rrig = 10 GeV) from the gluons (3.0%) principally, 
less from quarks (0.7%), independently of whether the Tevatron jet data are included in 
the fit. Since the jet data at large transverse energy Ej- are known to probe the behavior 
of g{x^Q'^) at large x, i.e., in the region where the depletion of the gluon's momentum is 
the strongest, we judge that inclusion of the jet data in the fit strengthens the constraining 
power of the fit. 
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FIG. 1: The gluon and singlet momentum distributions xg{x,Q^) and xT,{x,Q'^) are displayed 
as functions of x at Q = 10 GeV with as{Mz) = 0.118 and rrig = 10 GeV. The curves show the 
conventional CTEQ6M fit (solid) and the LG solutions with (short-dashed) and without (long- 
dashed) the Tevatron jet data included in the data set. 

C. Gluino contributions to hard scattering 

Once light superpartners are introduced as degrees of freedom, we must consider their 
impact on all hard scattering processes. Their effects can be felt both at tree level and 
in virtual-loop diagrams. At leading order in perturbation theory, we may consider hard 
subprocesses initiated by light gluinos or light bottom squarks that are constituents of the 
initial hadrons, as well as subprocesses in which gluinos or bottom squarks are emitted in 
the final state. We evaluate SUSY contributions to the hard matrix elements at leading 
order only for the same reasons that justify the omission of NLO SUSY contributions to the 
splitting kernels in Sec. IIIBI 

The CTEQ6 fit is performed to data from lepton-nucleon deep-inelastic scattering (DIS), 
vector boson production (VBP), and hadronic jet production at the Tevatron. In DIS and 
VBP, the lowest-order contribution from gluinos is 'y* + g ^ q + g at order 0{as). This 
subprocess proceeds via squark exchange, and its contribution can be neglected as being 
much smaller than the Born- level QCD subprocesses (which contribute at order 0{a^)). The 
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FIG. 2: The gluon, singlet, and gluino momentum distributions xg{x,Q'^), xT,{x,Q'^), and 
xg{x, Q"^) are displayed as functions of x at Q = 100 GeV with as{Mz) = 0.118 and rrig = 10 GeV. 
The curves show the conventional CTEQ6M fit (solid) and the LG solutions with (short-dashed) 
and without (long-dashed) the Tevatron jet data included in the data set. 

bottom squarks can contribute at C(a°) through the subprocesses 6 + 7* ^ 6 in DIS and 
6 + 6* — > 7* in the Drell-Yan process. However, these contributions appear in a combination 
with a small bottom squark PDF 6(x, Q^) and, therefore, are also negligible. We conclude 
that Born-level SUSY contributions appear only in the Tevatron jet production data, while, 
to the assumed level of accuracy, the hard matrix elements in DIS and VBP remain the 
same as in the SM case. 

We now consider the influence that gluino subprocesses may have on the rate for jet 
production at large values of transverse energy Et- Gluinos are color-octet fermions and, 
produced in the the final state, they materialize as jets. Since the gluino parton density is 
relatively large only at small values of fractional momentum x and, as illustrated in Figs. ^ 
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and 121 is small even there when compared with the gluon and light-quark densities, we are 
justified in neglecting the contributions to the large Et jet rate from subprocesses initiated by 
two gluinos. An example is g+g —* g+g. However, in the interest of completeness, we include 
two subprocesses initiated by one gluino: g + g ^ g + g, and q + g ^ q + g. Subprocesses 
initiated by gluons and/or light quarks can be important. We include g+g — ^ g + g via either 
a direct channel gluon or a cross channel gluino, and q + q ^ 9 + 9 via a direct channel gluon. 
We can ignore the t-channel exchange diagrams that contribute to q + q^g + g. With the 
possible exception of the bottom squark, the masses of most squarks are so large that the 
relevant t-channel amplitudes are negligible. In the case of bottom squark exchange, the 
two initial-state partons would be bottom quarks, with small parton densities. For similar 
reasons, we may also ignore subprocesses such as q + g -^ g + q. 

At the Tevatron pp collider, the qq partonic luminosity is relatively large in the region of 
large Ex, and one might expect naively that the subprocess q + q -^ g + g would increase the 
jet rate significantly. However, just as for the direct channel gluon subprocess q + q —* q' + q' 
in standard QCD, the squared matrix element for q + q—>-g + g is relatively small. 

In our treatment of jet production, we compute the matrix elements for the SUSY-QCD 
subprocesses at leading order. Working at large Et ^ rrig, we neglect the gluino mass in 
the calculation. We include these matrix elements as additional contributions to the jet rate 
in the fitting program, adding them to those of the NLO standard model QCD processes 
[0(a^)] to obtain constraints from the inclusive jet data. 

As indicated in the previous subsection, inclusion of a light gluino in the PDF set removes 
momentum from the gluon PDF at large x, tending to depress the contribution from SM 
processes to the jet rate at large Et- However, as we show, the effect is compensated 
partially by a larger value of as{Mz), by slower evolution of a.^ that makes as{ET > Mz) 
larger than in the standard model, and by contributions to the jet rate from production of 
SUSY particles in the final state. 

III. PRESENTATION OF THE GLOBAL FITS 

Our global fits are made to the complete set of data used in the CTEQ6 analysis, for 
several fixed values of the gluino mass. At this stage of the analysis, we do not impose a 
value of as{Mz)i preferring to determine a range of values directly from the global analysis 
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of hadronic scattering data. We perform fits to the hadronic data with as{Mz) set equal 
to one of several selected values in the range 0.110 < as{Mz) < 0.150. As discussed in 
Sec. I11C| jet production at the Tevatron is the only process among those we include that 
acquires Born-level contributions from the light gluinos. Most of the figures presented are 
for fits to the full set of data. To gauge the effect of the jet data, we also show results of 
additional fits that do not include these data. The numbers of experimental points in the 
fits with (without) the Tevatron jet data are 1811 (1688). 

A. Contour plots of A^^ vs a^ and rrig 

The principal result of our analysis is shown in Fig. El It maps the region of acceptable 
values of x^ iii the plane of as{Mz) and rrig. The contour plot shows the difference Ax^ = 
X^(as, rrig) — X^lcTE06M' between the value of x^ obtained in our LG fit and the standard 
model result equivalent to the CTEQ6M fit. The point in the plane corresponding to the 
CTEQ6M fit {as{Mz) = 0.118 and rrig -^ oo) is marked by the arrow. ^ The variation of x^ 
in the neighborhood of the minimum is used to estimate limits of uncertainty. 

An overall tolerance parameter T and a condition Ax^ < T^ are used in the CTEQ6 
analysis to characterize the acceptable neighborhood around the global minimum of x^ in the 
parton parameter space. The quantitative estimate T = 10 is obtained from a combination 
of the constraints placed on acceptable fits by each individual experiment included in the 

fitQ.' 

According to the tolerance on Ax^ of the CTEQ6 analysis, a fit is strongly disfavored if 
Ax^ > 100. The isoline corresponding to Ax^ = 100 is shown in Fig.Elby the solid line. The 
acceptable fits lie inside a trough that extends from large gluino masses and as{Mz) = 0.118 
down to rrig ^ 0.8 GeV and right to as{Mz) = 0.145. An even narrower area corresponds 
to fits with x^ close to those in the CTEQ6M fit. We note that x^ is better than in the 



^ The best-fit value as{Mz) = 0.117 in the CTEQ6 fit is slightly below the world-average value as{Mz) = 
0.118 assumed in the CTEQ6M PDF set. 

^ The tolerance T = 10 is estimated from the degree of consistency between the various data sets in the 
global fit. It includes effects due to experimental uncertainty and uncertainties that are of theoretical 
or phenomenological origin. It is an oversimplification to represent all uncertainties of PDFs and their 
physical predictions by a single number T. However, given the complexity of the problem, it is unrealistic 
to be more precise at this stage. The criterion T = 10 must be used with awareness of its limitations. 
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CTEQ6M fit in a small area in which rrig ^ 20 GeV and as{Mz) ^ 0.125, with the minimuni 
Ax^ ~ —25 at rrig = 8 GeV and as{Mz) = 0.130. This negative excursion in Ax^ is smaller 
than the tolerance T^ and should therefore not be interpreted as evidence for a light gluino. 

A substantial region of as{Mz) and rrig is excluded by the criterion Ax^ < 100. For 
as{Mz) = 0.118, gluinos lighter than 12 GeV are disfavored. However, the lower bound on 
rUg is relaxed to less than 10 GeV if as{Mz) is increased above 0.120. 

In Fig. El the positions are marked o 
PDF analyses with a light gluino 



the points {as{Mz), rrig} of the best fits in earlier 
^ Most of these earlier solutions are 



ed oi tne pomts 

HHQQ. 



the fits corresponding to rrig 



excluded by the present data set, with the exception o: 
GeV and large a,(M^) = 0.124, 0.129, and 0.134 lialnl- 

Another perspective is provided by a plot of x^ vs. rUg for several fixed values of as{Mz), 
shown in Fig. EJ The dependence on rrig is observed to be quasi-parabolic, with a shift of 
the minimum of x^ to lower rUg as as{Mz) increases. When as{Mz) is close to the current 
world-average value of 0.118, the fit prefers a heavy gluino, or no gluino at all. Very light 
gluinos are strongly disfavored, and the bound rrig > 12 GeV is obtained for Ax"^ < 100. 
For as = 0.122 and 0.124, the corresponding bounds are rrig > 5 GeV and rrig > 3 GeV, 
respectively. 

Conversely, for a very large as{Mz) (> 0.127), the pattern is reversed, and lighter gluinos 
[rrig < 50 GeV) are preferred. In the transition region of as{Mz) about 0.127, both very 
light and very heavy gluinos are disfavored, and a gluino mass in the range 10 to 20 GeV 
yields a slightly better x"^ than in the CTEQ6M fit. For a very small gluino mass of 1.3 to 
5 GeV, the minimum in x"^ is achieved for as{Mz) about 0.135, while even larger values of 
as{Mz) are disfavored as well (c.f. the curves for a^ = 0.140 and 0.145). 

The behavior of x^ in Figs. El and |3 exhibits irregular structure when the gluino mass lies 
in the range 50 to ~ 200 GeV. ^ These gluino masses lie beyond the range of sensitivity of 
the data sets in the fit, with the exception of the Tevatron jet data at jet transverse energies 
Et > 2mg. Variations in x^ at high gluino masses are caused by an interplay between the 
gluino mass and individual CDF and D0 jet data points. Contributions from gluinos with 
masses of 100 — 140 GeV improve the description of very high-ii^T- jet events, leading to a 



° The points corresponding to fits with gluino mass rrig < 0.7 GeV in Refs. |13jllallS| are off scale and are 

not shown. 
^ The irregularities of the contours are smoothed in Fig. |21 
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FIG. 3: A contour plot of Ax^ = x^ {cis{Mz),mg) — X^lcTE06M ^® ^ function of the strong cou- 
pling as{Mz) and gluino mass nig. The values of Ax^ are shown by labels on the corresponding 
isolines. The shaded region is exclude d by the CTEQ6 tolerance criterion. The points correspond- 



ing to the earlier PDF fits with a LG 



le 



13, 



IS 



19| are denoted by the symbols described in the 



legend. The solid line marks the Ax^ = 100 isoline. The dashed completion of this line at the 
bottom of the contour plot corresponds to fits done with rrig = Qo < nric = 1.3 GeV. 



dip in x^ in this region. Further discussion is found in Sec. 1111 C II 

In summary, as as{Mz) increases, the window of allowed gluino masses shifts from high 
to low values. Very large values of as(M^)(^ 0.148) can be ruled out for any gluino mass. 
Gluino masses between 10 and 20 GeV are allowed, as long as as{Mz) is not smaller than 
~0.118. 
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FIG. 4: Dependence of x^ on rrig for several fixed values of as{Mz) (denoted by labels on each 
corresponding line). 

B. Exploration of the light gluino fits 

The contour plot in Fig. El indicates that excellent fits to the global data can be obtained 
with a gluino mass below the weak scale, rrig < 100 GeV, provided that as(M^) is allowed 
to increase above the nominal value as{Mz) = 0.118. It is instructive to examine the 
compensating effects of as{Mz) > 0.118 and finite gluino mass on the parton distribution 
functions themselves. 

In Fig. 13 the ratio shown as a dashed line provides a comparison of the gluon distribution 
g{x,Q'^) at Q = 15 GeV and gluino mass nig = 15 GeV with g{x,Q^) in the CTEQ6M fit 
(without gluinos). The strong coupling as at the scale Mz is chosen to be the same as in the 
CTEQ6M fit, as{Mz) = 0.118. The ratio shows that, as the gluino mass is decreased below 
the weak scale, g{x, Q^) is depleted at large x and increased at small x. This softening of the 
gluon distribution follows from the slower evolution of ^^(Q), as well as from the presence 
of the additional coupling g ^ gg. 
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FIG. 5: Illustration of the compensation in the gluon density that arises between an increase in 
as{Mz) and a finite gluino mass. Shown are the ratios of the gluon density in the LG solution 
divided by the CTEQ standard model best fit (rrig = cxd and as{Mz) = 0.118) at Q = 15 GeV. 
The gluino mass is rrig = 15 GeV or rrig = oo, as specified in the figure. The CTEQ6 uncertainty 
band is shown by the dotted lines. 

For the same as{Mz), the magnitude of as{Q) at scales Q < Mz is smaller in the LG case 
than in the SM case (cf. Fig.ll2b). Correspondingly, PDF evolution is slower in the LG case. 
To some degree, the effects of the slower backward evolution can be compensated by selection 
of a larger value of as{Mz)- In some range of rrig and as{Mz), the effects of a smaller light 
gluino mass can be offset by a larger value of as{Mz)- This statement is illustrated by the 
dot-dashed curve in Fig. El The dot-dashed curve shows the ratio of the gluon density for 



rri!^ 



[m-n 



cxD and as{Mz) increased arbitrarily to 0.127, divided by the SM CTEQ6M result 
: cxD and as{Mz) = 0.118). The comparison shows that, when as{Mz) is increased, 
g{x, Q"^) is enhanced at large x and depleted at small x. The solid line in Fig. El shows that, 
by lowering rrtg below 100 GeV for a fixed ^^(M^), we can approximately cancel the effect 
of increasing as{Mz) at a fixed rrig. The solid line lies within the band of uncertainty of the 
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CTEQ6 gluon density, indicative of a fit of good quality. The cancellation breaks down at 
very large as{Mz)- 

A similar cancellation between the effect of a small rrig and increased as{Mz) is apparent 
in the singlet quark PDF. Consequently, a region exists at small rrig and large as{Mz) where 
the resulting PDFs remain close to those in the CTEQ6M fit. If as{Mz) is allowed to float 
freely in the fit, one can obtain PDFs similar to the CTEQ6M PDFs for all values of rrig 
above 0.8 GeV. For rrig ^ 150 GeV, the PDFs are practically the same as in the CTEQ6M fit, 
indicating that the current inclusive hadronic data are not sensitive to such heavy particles. 

C. Impact of various data sets 

To appreciate which data are the most restrictive in our fits, we examine the roles played 
in the fit by the hadronic jet data and other experiments. 

1. Tevatron jet data 

The Tevatron jet data places important constraints on rrig. In the absence of the jet data, 
the lower limit on rrig is weaker, with rrig ^ 5 GeV at as{Mz) = 0.118 if the jet data are 
omitted, but rrig > 12 GeV if the jet data are included. 

Comparisons between theor y an d the inclusive jet data from the CDF Collaboration 



and the D0 Collaboration 



34 



35| are shown in Fig. Inland Figs. [7|-|ni The results are from 
a series of fits for fixed as{Mz) = 0.118, 0.122, and 0.124. These three values of as{Mz) 
represent roughly the world-average central value and the values that are approximately 
one and two standard deviations larger. The CDF data are rescaled by 1.06, to account 
for differences in the measured luminosity used by the CDF and D0 Collaborations for the 
run-I data sample. The SM CTEQ6M fit with as{Mz) = 0.118 provides a good description 
of the data. For as{Mz) = 0.118 and rrig = 15 — 25 GeV, the theoretical cross sections 
are below the data due to the depletion of the gluon PDF at large momentum fractions x 
(Fig. EI). The gluon depletion can be counterbalanced in a wide range of rrig by a larger value 
of as{Mz)- The best fits correspond to combinations of a^ and rrig near the bottom of the 
trough in Ax^ in Fig. El For gluinos in the range 10-25 GeV, an acceptable fit is possible if 
as{Mz) is increased to about 0.124. 
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FIG. 6: Comparison of theoretical predictions with the CDF inclusive jet data. The data points 
show the ratio (Data-CTEQ6M)/CTEQ6M, where Data denotes the CDF measurements, and 
CTEQ6M is the SM prediction based on the CTEQ6M PDFs. The curves show the ratio (Theory- 
CTEQ6M)/CTEQ6M, where the Theory curves are the calculations based on the LG PDFs, for 
as{Mz) = 0.118, 0.122, and 0.124, and gluino masses m-g = {15,25,100} GeV. The solid curves 
correspond to the SM fits (effectively nig = oo) for the indicated choices of Qs(M^). The horizontal 
scale shows the transverse energy Et of the jet in GeV units. 

Contributions from gluinos increase the jet cross sections at Et > 2mg. A new channel 
for hard scattering is opened, and the evolution of as{Q) is slower. For as{Mz) = 0.118, a 
heavy gluino in the range 100 — 140 GeV improves agreement of theory with the Tevatron jet 
data in the high-ET tail by augmenting the rate of the tightly constrained standard model 
contributions. Better agreement for rrig = 100 GeV (dashed line) is visible in the high-ET 
region in Figs. Efa) and [7| Below the gluino threshold, the theory prediction (derived from 
the fit to the data insensitive to gluino contributions) is identical to the CTEQ6M fit. While 
X^ for the D0 data set is visibly improved (cf. Fig.llO|). the reduction of the overall x^ by 20 
units is not statistically significant. It will be interesting to see whether the trend in favor 
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FIG. 7: Comparison of theoretical predictions with the D0 jet data for as{Mz) = 0.118. Shown 
are (Data-CTEQ6M)/CTEQ6M and (Theory-CTEQ6M)/CTEQ6M, where Data are the D0 mea- 
surements, and CTEQ6M and Theory are defined as in Fig. |21 The data are sorted according to 
the ranges of the jet pseudorapidity ry. Curves are plotted for rrig = {15,25, 100} GeV. 

of contributions from gluinos (or other new color- charged fermions) with masses around 
100 GeV is preserved in the jet data from run-II at the Tevatron. 
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FIG. 8: Same as in Fig. for as{Mz) = 0.122. The solid curves correspond to the SM fit for 
as{Mz) =0.122. 

2. Plots of the data sets vs x^ 

To investigate further the influence of various sets of data, we display the ratios 
X^/XcTEQ6M fo^ individual experiments. Results for as{Mz) = 0.118 are shown in Fig. ITUl 
and those for as{Mz) = 0.135 in Fig. [TTl The choice of the extreme value as{Mz) = 0.135 
is made to accentuate the effects we want to demonstrate. In Fig. EH we observe that, in 
addition to the jet data (sets (n) and (o)), the DIS data from the HI Collaboration (sets 
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FIG. 9: Same as in Fig. for as{Mz) = 0.124. The solid curves correspond to the SM fit with 
as{Mz) = 0.124. 

(c) and (d)) and the CCFR F2 measurement (set (i)) tend to drive the gluino mass to large 
values when as{Mz) = 0.118. On the other hand, in Fig. [TT]we see that several sets of data 
are accommodated better with a light gluino if as{Mz) is large. 

If as{Mz) is chosen between 0.118 and 0.135, the behaviors of the values of x^ for indi- 
vidual experiments follow a mixture of the patterns shown in Figs. EH and ^B Several data 
sets disallow very small and very heavy gluino masses, while gluinos in the intermediate 
mass range are accommodated well by the fit. 
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FIG. 10: Ratios X^/Xcteq6M ^^ individual experiments and and as{Mz) = 0.118. The curves 
correspond to the following data sets: (a) BCDMS F| [3^; (b) BCDMS F^ ^]; (c) HI F^ 



(1996/97) 
NMC Fi/FF 
production 
production 



39|; (d) HI Ff (1998/99) |40|; (e) ZEUS F| (1996/97) |4l|; (f) NMC Ff '42\; (g 



4|; (h) NMC Fi/F^ |4^; (i) CCFR Fa |4j|; (j) CCFR xFg ^A% (k) E605 muon pair 

45|l; (1) CDF lepton asymmetry [4a|; (m) E866 muon pair production .47.] : (n) D0 jet 

35l |: and (o) CDF jet production 



D. Section summary and momentum fractions 

If as{Mz) is allowed to vary freely, reasonable fits to the global data set are possible for 
essentially any gluino mass above ~ 1 GeV. However, if as{Mz) is constrained from other 
sources, say, r decay or direct measurements at Mz, then a global fit to scattering data im- 



poses good constraints on rrig. This situation is reminiscent of the strong corre 
the gluon PDF and a^ observed in previous analyses of parton densities 



ation between 



49|. Similarly, 



it is not surprising that constraints on the gluino mass are coupled to our knowledge of the 



22 



1.5 



0.5 
1.5 



0.5 



X^irn^MXcimm foras[Mz]=0.135 




mg,GeV 



mg, GeV 



FIG. 11: Same as in Fig-CUl for as{Mz) = 0.135. 

gluon PDF (constrained by the hadronic jet data) and as{Mz)- 

The principal uncertainties on our quoted bounds on the gluino mass arise from neglect of 
NLO supersymmetric contributions to the PDF evolution (affecting the PDFs at a percent 
level), neglect of NLO SUSY-QCD corrections to jet production, and the hmited precision 
of the criterion Ax^ < 100 for the selection of acceptable fits. The lower limit on the gluino 
mass can be relaxed if the NLO virtual-loop SUSY-QCD corrections enhance the rate of 
the standard model subprocesses in the Tevatron jet production. These uncertainties can 
be reduced in future analyses. 

We conclude this section with Table HI in which we show the fraction of the proton's 
momentum carried by its constituents, in both the standard model CTEQ6M fit and in the 
LG fit with nig = 15 GeV and a^^Mz) = 0.122. 
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Parton type 


CTEQ6M 


nig = 15 GeV, 
as{Mz) = 0.122 


u + u 


25.2 


25.4 


d + d 


15.4 


15.5 


s + s 


5.3 


5.1 


c + c 


3.9 


3.6 


b + b 


2.5 


2.3 


S 


52.3 


51.9 


9 


47.5 


44.7 


9 





3.2 


Total: 


100 


100 



TABLE I: Fractions of the proton's momentum carried by different parton species at Q = 100 GeV 
in the CTEQ6M fit {rrig -^ oo and as{Mz) = 0.118) and in the LG fit with mass nig = 15 GeV 
and as{Mz) = 0.122. 

IV. SUMMARY 



Our new analysis of the com pat ibility of a light gluino with inclusive scattering data goes 
beyond earlier studies jla, IhI, UM, u^ in a number of aspects. First, the current data are 
strikingly more extensive than available ten years ago. They cover both small-x and large-Q 
regions, come from a variety of experiments, and are characterized by high precision. The 
primary effects of a gluino in the global analysis are changes in the evolution of the strong 
coupling strength and changes in the evolution of the parton distributions. It is easier to 
observe these changes in a data sample with large lever arms in Q and x. 

Second, in contrast to previous studies, our fit includes the complete set of data from 
the CTEQ global analysis, including the Tevatron jet production data. The major role of 
the jet data is to constrain the gluon density at large values of fractional momentum x. 
The behavior of the gluon density at large x is affected strongly by the presence of the 
gluinos in the mix. Because they are sensitive to gluons at large x, the jet data enhance the 
discriminating power of the global fit. When the gluino mass changes, large variations in 
X^ are observed, an influence that can be used to constrain the gluino parameter space. In 



24 



view of the strong correlations between the gluino mass, as{Mz), and the gluon distribution, 
the constraints can be determined only after a consistent implementation of SUSY effects 
throughout all stages of the analysis. 

The third new component in our study is a method [7] for quantitative interpretation 
of uncertainties in parton distributions. With the help of this method, constraints on the 
acceptable gluino mass can be imposed on the basis of the values of x^ obtained in the fits. 
The main result of the paper is presented in Fig. El It shows the region of the gluino masses 
rrig and QCD coupling strengths as{Mz) allowed by the present data. The standard model 
fit prefers as{Mz) ~ 0.118. Gluinos with mass of a few GeV can be accommodated only if 
as{Mz) is increased. For example, gluinos with mass below 1 GeV are admissible only if 
as{Mz) is about 0.130 or larger. If one takes into account the effects of a light gluino on 
the measurement and running of a^, it is hard (if at all possible) to reconcile such a large 
value of as{Mz) with both low- and high-energy electroweak data. 

On the other hand, a possibility remains open for the existence of gluinos with mass 
between 10 to 20 GeV, with a moderately increased as {as{Mz) > 0.119). This possibility 
is even slightly favored by the x^ of the global fit. A model with light gluinos and bottom 
squarks 6] is not incompatible with the results of our PDF analysis. Tighter constraints 
can be obtained in the near future, when new data from HERA and the Tevatron become 
available. In particular, it will be intriguing to see whether higher statistics jet data at 
larger values of Et show indications of physics beyond the standard model. The constraints 
we obtain depend on the value of as{Mz). Uncertainties in the value of as{Mz) could 
be reduced through a consistent determination of as{Mz) from the CERN LEP data in a 
SUSY-QCD analysis, in which the effects of superpartners are included in the data analyses 
and Monte-Carlo simulations. 

Implementation of the gluino in our study relies only on the knowledge of its strong 
interactions, which are determined uniquely by supersymmetry. We consider only theoret- 
ically clean one-scale inclusive observables. In this sense, our constraint rrig > 12 GeV for 
as{Mz) = 0.118 should be compared to the constraint rrig > 6.3 GeV for the same value 
of as from the Z-boson width measurement ^]. Tighter constraints on rrig were quoted 
by the searches for traces of gluino hadronization J5l|, |52| and a study of jet shapes J53[. 
Although important, these constraints are less general, since they involve assumptions about 
the gluino lifetime or deal with several momentum scales in the jet shape observables. Our 
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TABLE II: as{Mz) derived by evolution from as(m^) = 0.323±0.030 [2J1 for several gluino masses 
and in the standard model. 



rrig, GeV 


as{Mz) 


10 


0.132 ±0.004 


25 


0.125 ±0.004 


50 


0.121 ±0.004 


90 


0.118 ±0.003 


SM 


0.118 ±0.003 



study demonstrates the potential of global analysis to independently constrain new physics 
from hadron collider data. 
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APPENDIX: STRONG COUPLING STRENGTH 



In this appendix we discuss the quantitative changes that occur in the evolution of the 
strong coupling strength as{Q) in the presence of a light gluino in the spectrum. We consider 
compatibility of r decay and LEP Z-pole measurements of as, which constrain as{Q) at low 
and large momentum scales, respectively. 

The measurement of a^ at LEP provides a constraint at scales of order of the Z-boson 
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FIG. 12: Dependence of as{Q) on the renormalization scale Q. The three data points shown are 



2J], and a<,(M2) = (0.118- 0.126) ±0.005 |27|. 



asirur) = 0.35 ±0.03 _54], as{mr) = 0.323 ±0.030 
In Fig. 112b . the average asirrir) = 0.337 of two experimental values at Q = rrir is evolved to higher 
energies. In Fig. 112b . the central value as{Mz) = 0.122 of the interval at Q = Mz is evolved to 
lower energies. The thick solid line represents the Standard Model (SM) evolution in the absence 
of SUSY effects. The dashed series of curves are generated for gluino masses Tiig = 5, 10, and 25 
GeV (shown by the labels on the corresponding curves). 

mass Mz, and the measurement of as in r-lepton decay provides a constraint at scales of 
order of the r-lepton mass 171^. If the gluino is substantially heavier than the r-lepton, 
its presence in the spectrum cannot affect the measurement of a^ in r decay. Therefore, 



as{Q) measured at the scale Q = m,- i s t 
quoted values are 0^(771,-) = 0.35 ± 0.03 



^e same as in the standard model. Its recently 
and asirrir) = 0.323 ± 0.030 |2J . On the other 
hand, the measurement of as{Q) at Q = Mz can be affected by light superpartners in the 
spectrum. The values obtained in standard model analyses may have to be revised. A recent 
estimate 271] of SUSY effects on direct measurements of as{Mz) finds values in the interval 
(0.118 — 0.126) ± 0.005, where the variation in the central value arises from the uncertainty 
in the magnitude of SUSY-QCD corrections to the Z-boson decay width from processes such 
as Z ^ bbg and Z -^ bbg. 

The r-decay and LEP values of as{Q) must be related by the renormalization group 
equation. Figure IT^ shows the evolution of as{Q) measured in r decay to the energy of order 
Mz for different choices of the gluino mass. According to Eqs. Q and Q, «s(Q) evolves 
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more slowly in the presence of light gluinos. Forward evolution of as {mr) in the LG case leads 
to a higher value at Q = Mz than in the standard model. Evolution of as{fnr) = 0.35 ±0.03 
results in as{Mz) = 0.120 ±0.003 in the standard model and 0.135 ±0.004 for a gluino mass 
of 10 GeV. Table im lists the values of as{Mz) obtained by evolution from the r decay value 
asirrir) = 0.323 ± 0.030. 

Alternatively, we may evolve as measured in Z-boson decay backward to energies of order 
rrir (Fig. IT^). The resulting asirrir) in the LG case is lower than in the standard model. 



If we use the central value as{Mz) = 0.122 from the interval 0.118 — 0.126 of Ref. J27| as 
our starting point, we obtain as{mr) = 0.367 in the standard model and 0.266 for rrig = 10 
GeV. 

For the quoted experimental and theoretical uncertainties, the measurements of a^ at rrir 
and Mz are incompatible at the one-standard deviation (1 a) level if the gluino mass is less 
than about 5.8 GeV. For such gluino masses, the 1 a interval of as{Mz) obtained by evolution 
from asifrir) is above the 1 a interval for the LEP measurement. On the other hand, the r 
decay and LEP data do agree at the 1 a level for a gluino heavier than 5.8 GeV, if as{Mz) is 
at the upper end of the theoretical uncertainty range (i.e., as{Mz) ~ 0.126 ± 0.005). Lower 
values of as{Mz) increase the lower bound on rrig, but gluino masses in the range 10 — 25 
GeV are possible within the uncertainties, as long as the central value of as{Mz) from LEP 
is not less than about 0.119. 

Bounds on the gluino mass obtained from the global PDF fit depend on the value of 
as{Mz) assumed in the fit. Gluino masses lighter than 6 GeV cannot agree simultaneously 
with the results of the global fit, r decay, and LEP Z pole measurements. Gluino masses 
as small as 10 GeV are consistent with values of as{Q) obtained from both r decay and 
LEP data, if as{Mz) is increased moderately compared to its SM world-average value of 
0.118 ±0.002. 
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